Calibration of Thermal Viscoelastic Material Models for the Dynamic Responses of PVB and SG Interlayer Materials

Laminated glass interlayer materials polyvinyl butyral (PVB) and SentryGlas® (SG, kuraray, Houstan, TX, USA) exhibit thermal viscoelastic behavior under dynamic tensile loading. Significant temperature and strain rate effects on the behavior of these interlayer materials pose a challenge for accurately modeling the dynamic response of laminated glass. Many researchers have simplified their approaches by modeling the response of the interlayer material using a bilinear approximation or established hyperelastic models. However, temperature and strain rate effects can be captured using the three-network viscoplastic (TNV) model. Therefore, the objective of this study is to calibrate material models for the thermal viscoelastic dynamic responses of PVB and SG interlayer materials. Uniaxial tensile tests were performed at strain rates of 2, 20, and 45 s−1 and temperatures of 0, 23, and 60 °C, and material models were calibrated using the experimental data. Finite element analysis using the calibrated material models successfully predicted the dynamic responses of PVB and SG under the experimental test conditions within a 10% error margin. This suggests that the calibrated models using the TNV model represent significant improvements over existing approaches to modeling the dynamic response of laminated glass. Similar procedures can be applied to other thermoplastics, laying the groundwork for establishing a standard calibration guide.


Introduction
The surge in terrorist attacks targeting densely populated urban areas in recent years has heightened concerns about building safety, underscoring the critical importance of designing building envelopes, which increasingly incorporate large glass curtain wall systems.In the event of failure in such systems, sharp glass fragments may become projectiles, posing a serious threat to occupants and potentially causing injuries or fatalities [1].Studies indicate that the majority of nonfatal injuries sustained during bomb blasts stem from flying fragments originating from architectural glazing [2,3].Laminated glass (LG) has been found to be effective at mitigating these risks by retaining glass fragments on a polymeric interlayer upon fracture [4][5][6].Additionally, after glass cracking, the ductile polymer interlayer deforms significantly as a continuous membrane providing post-cracking energy absorption [7].The shear behavior of the interlayer material affects the deformability and load-bearing capacity of laminated glass plates [8], as well as their post-breakage behavior [9][10][11].In other words, the polymeric interlayer governs the general behavior of laminated glass, highlighting its importance for the adequate design and calculation of laminated glass structural elements.Interlayer materials are viscoelastic, meaning their stiffness and mechanical response depend on load duration and working temperature [12][13][14].Therefore, the mechanical properties of laminated glass also depend on at least these two factors [8,15,16].
Larcher et al. [17] state that the failure process of a laminated glass sheet can be divided into five phases: (1) elastic behavior of glass plies, (2) outer glass ply cracks while the interlayer is not damaged, (3) inner glass ply fails and the interlayer behaves elastically, (4) interlayer deforms plastically and splinters are kept together by the interlayer, and (5) the interlayer fails at ultimate strength or is cut by the glass shards.Extensive research has been conducted on the elastic behavior of laminated glass, making modeling phase (1) with numerical models straightforward.However, phases (3) through (5), largely dominated by the mechanical properties of the viscoelastic polymer interlayer, pose greater challenges.Thus, calibrating material models to accurately predict the dynamic responses of polymeric interlayer materials is critical for the analysis and design of laminated safety glass.
Kott and Vogel define the failure modes of laminated glass panels to specify their remaining structural capacity [18].These preconditions, along with four-point bending tests and previous impact tests, form the foundation for developing a new mechanical model.Various geometric parameters are evaluated for panes supported on two sides, considering both limit states: load-carrying capacity (LCC) and residual strength capacity (RSC).The types of fractures and the corresponding yield line mechanisms are also analyzed.This comprehensive approach enables the improvement of RSC and provides an analytical verification.
In the literature, various approaches have been explored for modeling laminated glass windows due to their inherent complexity as composite materials composed of constituent materials with vastly different mechanical properties [6,39].Among the finite element models utilized for laminated glass, researchers often combine shell elements with membrane elements [17,[40][41][42][43]. Glass is typically modeled using shell elements as linear elastic materials with failure at a critical strain limit [17,[40][41][42][43] or Rankine stress limit [42], or with piecewise linear plasticity models to allow small amounts of plastic strain [40,43].However, modeling of the interlayer material varies.Some researchers have approximated the polyvinyl butyral (PVB) interlayer as linear elastic, neglecting large nonlinear strain behavior and viscoelasticity [43].Others have employed widely used hyperelastic laws available in LS-DYNA, such as the Blatz-Ko, Mooney-Rivlin, and Ogden materials, to model the PVB interlayer [40,41,44].While hyperelastic material models capture the nonlinear elastic behavior of materials based on their stress-strain responses, they do not inherently account for strain rate or temperature dependencies and do not accurately capture the stiff initial linear elastic region of polymers under high strain rates.
Larcher et al. [17] noted that the behavior of PVB under higher strain rates is more elastic with hardening than viscoelastic, and accordingly modeled PVB interlayer with a bilinear elastic-plastic material law, as shown in Figure 1.Lusk et al. [45] modeled laminated glass windows under blast loading using the LS-DYNA material model 32 "MAT-LAMINATED_GLASS", which treats the response of the interlayer as bilinear with input data of Young's modulus, Poisson's ratio, yield stress, and plastic hardening modulus.Many other researchers have employed this elastic-plastic bilinear approach to modeling the response of laminated glass interlayer materials [37,[46][47][48][49].However, while the elastoplastic material model approach approximates the stiff initial linear region at high strain rates, it does not account for the significant strain rate and temperature dependencies in the mechanical responses of interlayer materials.Furthermore, it does not accurately capture more complex material responses, such as that of the interlayer material SentryGlas ® (SG).
modeling the response of laminated glass interlayer materials [37,[46][47][48][49].However, while the elasto-plastic material model approach approximates the stiff initial linear region at high strain rates, it does not account for the significant strain rate and temperature dependencies in the mechanical responses of interlayer materials.Furthermore, it does not accurately capture more complex material responses, such as that of the interlayer material SentryGlas ® (SG).[36], Iwasaki [50], Bennison [22]).
Analytical approaches have been made to develop master curves for the mechanical characterization of viscoelastic materials that incorporate time and temperature variables.Hána et al. formulated master curves for two types of interlayers, such as PVB-based TROSIFOL ® BG R20 and EVA-based EVALAM ® , using the results of static and dynamic single-lap shear tests to identify parameters for a generalized Maxwell chain model with an incorporated temperature shift factor to describe their time-and temperaturedependent relaxation functions [51].Similarly, Centelles et al. performed relaxation tests on seven polymeric films used as laminated glass interlayer materials and developed master relaxation curves using a generalized Maxwell model and the time-temperature superposition principle [52].Liu et al. [29] performed tensile tests on PVB material under a range of strain rates and proposed sets of constitutive equations using the ZWT model.The parameters of the constitutive equations were fitted to each individual strain rate and tested without considering temperature effects.Chen et al. performed uniaxial tensile tests on the ionomer interlayer (SG) at a wide range of strain rates and temperatures and proposed dynamic constitutive equations for different temperatures based on a modified G'Sell model [24].The results of these analytical approaches did not yield material models readily available for implementation into finite element analysis and were not validated by numerical simulations.
The mechanical responses of thermoplastics, such as PVB and SG, are typically well captured with a three-network frame [53].Two efficient candidates are the three-network (TN) model and the three-network viscoplastic (TNV) model [54].Both the TN and TNV models are powerful models capable of capturing strain rate and temperature effects.The three-network (TN) model is specifically developed for thermoplastic materials, and as such, the individual networks implemented are immutable.Shahin et al. [55] calibrated the TN model to successfully model the relaxation response of high-density polyethylene (HDPE) flanges with an error less than 10%.Similar model calibrations for HDPE that are solely suitable for constant loads, e.g., creep relaxation, were proposed by Lai et al. [56] Figure 1.Elastic-plastic material law used by Larcher et al. [8] compared with experimental test data of PVB at different strain rates (Morrison [36], Iwasaki [50], Bennison [22]).
Analytical approaches have been made to develop master curves for the mechanical characterization of viscoelastic materials that incorporate time and temperature variables.Hána et al. formulated master curves for two types of interlayers, such as PVB-based TROSIFOL ® BG R20 and EVA-based EVALAM ® , using the results of static and dynamic single-lap shear tests to identify parameters for a generalized Maxwell chain model with an incorporated temperature shift factor to describe their time-and temperature-dependent relaxation functions [51].Similarly, Centelles et al. performed relaxation tests on seven polymeric films used as laminated glass interlayer materials and developed master relaxation curves using a generalized Maxwell model and the time-temperature superposition principle [52].Liu et al. [29] performed tensile tests on PVB material under a range of strain rates and proposed sets of constitutive equations using the ZWT model.The parameters of the constitutive equations were fitted to each individual strain rate and tested without considering temperature effects.Chen et al. performed uniaxial tensile tests on the ionomer interlayer (SG) at a wide range of strain rates and temperatures and proposed dynamic constitutive equations for different temperatures based on a modified G'Sell model [24].The results of these analytical approaches did not yield material models readily available for implementation into finite element analysis and were not validated by numerical simulations.
The mechanical responses of thermoplastics, such as PVB and SG, are typically well captured with a three-network frame [53].Two efficient candidates are the three-network (TN) model and the three-network viscoplastic (TNV) model [54].Both the TN and TNV models are powerful models capable of capturing strain rate and temperature effects.The three-network (TN) model is specifically developed for thermoplastic materials, and as such, the individual networks implemented are immutable.Shahin et al. [55] calibrated the TN model to successfully model the relaxation response of high-density polyethylene (HDPE) flanges with an error less than 10%.Similar model calibrations for HDPE that are solely suitable for constant loads, e.g., creep relaxation, were proposed by Lai et al. [56] and Elluech et al. [35,57].The three-network viscoplastic (TNV) model is more versatile and is capable of capturing the time dependence, the pressure dependence of plastic flow, the pressure-dependent bulk modulus, the volumetric plastic flow, and anisotropy.The individual networks making up the TNV model can be selected and optimized to fit a wide range of material behaviors.Recently, Almomani et al. [53] calibrated the TNV model to predict the thermal-viscoelastic behaviors of HDPE under monotonic and cyclic loads.The model was implemented into FEA analysis and was successfully able to predict the HDPE behavior under static, quasi-static, and dynamic loads and at temperatures from 23 • C to 95 • C.
Based on the aforementioned studies, it can be concluded that the TNV model shows promise for capturing the versatile, nonlinear strain rate and temperature-dependent behaviors of laminated glass interlayers.However, to date, no studies have been found that investigate its calibration on PVB, SG, or any laminated glass interlayer materials for the monotonic tensile response at various strain rates and temperatures.Therefore, the objective of this study is to calibrate the TNV model to capture the thermo-viscoelastic dynamic responses of PVB and SG interlayer materials under a wide range of strain rates and temperatures.The calibrated material models will then be validated by finite element simulations.

Experimental Program
In this section, the interlayer materials tested are described in detail.Test specimen preparations are provided.An overview of two tensile test methodologies, pneumatic testing and drop-weight testing, is presented.Data acquisition and processing are discussed.Finally, the experimental test matrix is detailed.

Material Description
Polyvinyl butyral (PVB), a solid thermoplastic resin invented in 1927, was the first material used as an interlayer for laminated glass, and it remains widely used today in various structural applications.PVB is produced through the reaction of polyvinyl alcohol with butyraldehyde.This random amorphous polymer consists of three monomers, each contributing specific properties: hydrophobic and elastic vinyl butyral (~80 wt%) provides processability, while hydrophilic vinyl alcohol (~18 wt%) and vinyl acetate (~2 wt%) ensure high adhesion to inorganic materials such as glass [58].
Several experiments have been conducted to understand the dynamic response of PVB.Nawar et al. showed that the mechanical response of PVB is highly strain-rate-dependent and that the quasi-static behavior resembles a hyperelastic material such as rubber, and at high strain rates, the behavior transitions to a viscoelastic bilinear response [37].Andreozzi et al. performed experiments with a rheometer and found that a PVB between 0 • C and 50 • C is in the glass transition region.Therefore, its dynamic modulus is highly temperaturedependent in the typical in situ temperature ranges.At high temperatures, the dynamic modulus is extremely low, and the behavior is hyperelastic.At low temperatures, PVB behaves like an elasto-plastic material and exhibits some brittleness [26].
SentryGlas ® (SG), sometimes referred to as an "ionoplast interlayer", is the only commercially available ionomer-based interlayer material for laminated glass.Developed by DuPont for use in structural safety glazing, SG offers much greater rigidity, higher strength, and better ductility compared to PVB [24].Since its market release in 1998, structural engineers have commonly used SG as a replacement for PVB to improve the mechanical response of laminated glass.
An ionomer is a polymer that typically has an ionic content of no more than 10 mole percent.The general chemical structure of an ionomer consists of a hydrocarbon containing neutralized pendant acid groups [59].SG is formed through a process in which a copolymerization of ethylene with methacrylate is cured with metal ions [60].Due to its semi-crystalline chemical structure and crosslinking, SG offers superior mechanical properties compared to traditional amorphous polymer interlayers.Furthermore, SG provides excellent adhesion to glass and metals because of its ionic interactions and chemical bonding.
Several studies have been performed to investigate the thermo-viscoelastic properties of SG interlayer materials.Zhang et al. [23] performed uniaxial tensile tests on SentryGlas Plus (SGP) interlayer material under a wide range of strain rates and found that SGP exhibited an elasto-plastic behavior at all strain rates, ranging from quasi-static rates to extremely high dynamic rates.Due to its semicrystalline structure, SG undergoes a yielding event after a stiff initial linear region and subsequent strain softening as the polymer chains begin to slide past each other within the amorphous regions and the crystalline domains continue to deform elastically.The strain rate effect predominantly affects the initial response, resulting in increased yield stress with increasing strain rate.The tests also found that SGP material became less ductile as the strain rate increased.Chen et al. [24] performed uniaxial tensile tests on SG interlayer material at a wide range of strain rates and temperatures and found that the initial modulus decreased with the increase in temperature; the yield stress decreased sharply with the increase in temperature, especially at high strain rates; and the yield strain showed limited dependence on the strain rate.The glass transition temperature of SG was approximately 50-55 • C [25].
Laminated glass manufacturers typically receive interlayer materials in sheets, with specific thicknesses of rolls.The sheets are then cut to size and placed in between the glass plies, and the composite window is processed in an autoclave under high pressures and temperatures.The curing process has been shown to have significant effects on the material responses of laminated glass interlayers.Accordingly, PVB and SG sheets that had already gone through the curing process were acquired.

Specimen Preparations
In this work, the sample geometry for the intermediate strain rate testing at a 2 s −1 strain rate, performed utilizing the pneumatic testing machine, was the ASTM D638-14 [61] standard Type IV specimen, as shown in Figure 2a.A modified ASTM D638-14 standard Type I geometry was employed for the high-strain-rate testing at 20 s −1 and 45 s −1 strain rates, which was performed utilizing a drop-weight testing machine.The end area was increased, as shown in Figure 2b, and aluminum tabs were adhered to the end tabs of the specimens to prevent tearing-out failure of the specimen under impact loading.Steel cutting dies were fabricated for both specimen geometries to ensure precise dimensions [38].For each specimen, the 1 inch central gauge length, L g , was marked with thin black lines/points using a permanent marker pen prior to testing to enable the strain to be monitored during the test using a high-speed camera.Steel digital calipers were used to measure the thickness and width of the test section at three locations, to an accuracy of 0.0005 in.chemical bonding.
Several studies have been performed to investigate the thermo-viscoelastic properties of SG interlayer materials.Zhang et al. [23] performed uniaxial tensile tests on SentryGlas Plus (SGP) interlayer material under a wide range of strain rates and found that SGP exhibited an elasto-plastic behavior at all strain rates, ranging from quasi-static rates to extremely high dynamic rates.Due to its semicrystalline structure, SG undergoes a yielding event after a stiff initial linear region and subsequent strain softening as the polymer chains begin to slide past each other within the amorphous regions and the crystalline domains continue to deform elastically.The strain rate effect predominantly affects the initial response, resulting in increased yield stress with increasing strain rate.The tests also found that SGP material became less ductile as the strain rate increased.Chen et al. [24] performed uniaxial tensile tests on SG interlayer material at a wide range of strain rates and temperatures and found that the initial modulus decreased with the increase in temperature; the yield stress decreased sharply with the increase in temperature, especially at high strain rates; and the yield strain showed limited dependence on the strain rate.The glass transition temperature of SG was approximately 50-55 °C [25].
Laminated glass manufacturers typically receive interlayer materials in sheets, with specific thicknesses of rolls.The sheets are then cut to size and placed in between the glass plies, and the composite window is processed in an autoclave under high pressures and temperatures.The curing process has been shown to have significant effects on the material responses of laminated glass interlayers.Accordingly, PVB and SG sheets that had already gone through the curing process were acquired.

Specimen Preparations
In this work, the sample geometry for the intermediate strain rate testing at a 2 s −1 strain rate, performed utilizing the pneumatic testing machine, was the ASTM D638-14 [61] standard Type IV specimen, as shown in Figure 2a.A modified ASTM D638-14 standard Type I geometry was employed for the high-strain-rate testing at 20 s −1 and 45 s −1 strain rates, which was performed utilizing a drop-weight testing machine.The end area was increased, as shown in Figure 2b, and aluminum tabs were adhered to the end tabs of the specimens to prevent tearing-out failure of the specimen under impact loading.Steel cutting dies were fabricated for both specimen geometries to ensure precise dimensions [38].For each specimen, the 1 inch central gauge length, Lg, was marked with thin black lines/points using a permanent marker pen prior to testing to enable the strain to be monitored during the test using a high-speed camera.Steel digital calipers were used to measure the thickness and width of the test section at three locations, to an accuracy of 0.0005 in.

Test Setups
For intermediate strain rate testing, a pneumatic testing machine was designed and fabricated to be able to perform temperature-controlled tensile tests at strain rates ranging from quasi-static rates to 5 s −1 , as shown in Figure 3a.The device was equipped with a load cell with a capacity of 250 pounds (1112 N) and was operated under displacement control by a LabView program.The elongation of the specimens throughout testing was tracked using an edgertronic SC1 model high-speed camera, capturing frames at a rate of 3000 frames per second.Select frames from a video of an intermediate strain rate test on PVB are shown in Figure 3b.A temperature enclosure was built around the testing frame to enable us to perform temperature-controlled testing.The high-strain-rate tensile testing in this research was performed using a m drop-weight apparatus; see Figure 4.This machine was originally invented for testing of robust materials like metals or composites.For the performance of high rate examinations on polymers, multiple novel components were designed, fab and implemented.An illustrative diagram of the machine, with the critical parts ide is presented in Figure 4.A moveable member with variable weight, shown in Fi was equipped with forked strikers, as shown in Figure 4, and installed onto the guid A fixed member with bore holes was fabricated to allow the strikers to pass throu impact the lower clamp, which was attached to the bottom tab of the test specimen, as to stop the moveable member after the specimen had failed.The upper tab of th imen was attached to a load cell with the upper clamp.A piezoelectric load cell (O dyne model LC213-500; Omegadyne, Nowalk, CT, USA; Edgertronic, Campbell, CA was used to record the load at a rate of 3000 data points per second.The elongation specimens throughout testing was tracked using an edgertronic SC1 model high camera, capturing frames at a rate of 3000 frames per second.The load from the lectric load-cell was recorded using a National Instruments USB-6351 data acquisiti tem.The high-strain-rate tensile testing in this research was performed using a modified drop-weight apparatus; see Figure 4.This machine was originally invented for impact testing of robust materials like metals or composites.For the performance of high-strainrate examinations on polymers, multiple novel components were designed, fabricated, and implemented.An illustrative diagram of the machine, with the critical parts identified, is presented in Figure 4.A moveable member with variable weight, shown in Figure 4, was equipped with forked strikers, as shown in Figure 4, and installed onto the guide rails.A fixed member with bore holes was fabricated to allow the strikers to pass through and impact the lower clamp, which was attached to the bottom tab of the test specimen, as well as to stop the moveable member after the specimen had failed.The upper tab of the specimen was attached to a load cell with the upper clamp.A piezoelectric load cell (Omegadyne model LC213-500; Omegadyne, Nowalk, CT, USA; Edgertronic, Campbell, CA, USA) was used to record the load at a rate of 3000 data points per second.The elongation of the specimens throughout testing was tracked using an edgertronic SC1 model high-speed camera, capturing frames at a rate of 3000 frames per second.The load from the piezoelectric load-cell was recorded using a National Instruments USB-6351 data acquisition system.
A temperature enclosure was built around the testing frame to be able to perform temperature-controlled testing.The temperature control system utilized was comprised of a Thermo Fisher Scientific™ Accel 500 LT cooling/heating chiller with a working temperature range of −25 • C to +80 • C with an accuracy of ±0.5 • C, an industrial-grade air ventilator, and insulated ducts.The temperature control system with the drop weight testing machine is shown in Figure 5.A temperature enclosure was built around the testing frame to be able to perform temperature-controlled testing.The temperature control system utilized was comprised of a Thermo Fisher Scientific™ Accel 500 LT cooling/heating chiller with a working temperature range of −25 °C to +80 °C with an accuracy of ±0.5 °C, an industrial-grade air ventilator, and insulated ducts.The temperature control system with the drop weight testing machine is shown in Figure 5.

Data Processing
The load-time history for each test, measured by the load cell, was recorded in Lab-View programs at 3000 samples per second.The engineering stress was calculated by dividing the load by the original cross-sectional area.Video of each test was captured by a high-speed camera at 3000 frames per second and later converted to files of sequential JPEGs and post-processed using proprietary software developed in-house.The engineering strain-time history of each test was produced by tracking the elongation of the central gauge length dot/line markings made on the specimens prior to testing.For each test, the  A temperature enclosure was built around the testing frame to be able to perform temperature-controlled testing.The temperature control system utilized was comprised of a Thermo Fisher Scientific™ Accel 500 LT cooling/heating chiller with a working temperature range of −25 °C to +80 °C with an accuracy of ±0.5 °C, an industrial-grade air ventilator, and insulated ducts.The temperature control system with the drop weight testing machine is shown in Figure 5.

Data Processing
The load-time history for each test, measured by the load cell, was recorded in Lab-View programs at 3000 samples per second.The engineering stress was calculated by dividing the load by the original cross-sectional area.Video of each test was captured by a high-speed camera at 3000 frames per second and later converted to files of sequential JPEGs and post-processed using proprietary software developed in-house.The engineering strain-time history of each test was produced by tracking the elongation of the central gauge length dot/line markings made on the specimens prior to testing.For each test, the

Data Processing
The load-time history for each test, measured by the load cell, was recorded in LabView programs at 3000 samples per second.The engineering stress was calculated by dividing the load by the original cross-sectional area.Video of each test was captured by a high-speed camera at 3000 frames per second and later converted to files of sequential JPEGs and postprocessed using proprietary software developed in-house.The engineering strain-time history of each test was produced by tracking the elongation of the central gauge length dot/line markings made on the specimens prior to testing.For each test, the failure stress (σ f ) and corresponding strain at failure (ε f ) were calculated.The initial Young's modulus, E 0 , was evaluated as the slope of the first portion of the linear response.The yield stress (σ y ) and corresponding yield strain (ε y ) were taken at the local maximum stress after the initial linear region.The tensile toughness, U, was calculated as the total area under the stress-strain curve.

Experimental Testing Matrix
PVB (RA41) from Eastman and SG5000 from Kuraray were tested at temperatures of 0, 23, and 60 • C and strain rates of 2, 20, and 45 s −1 .The experimental testing matrix is given below in Table 1.

Experimental Results
In this section, the engineering stress-strain results for the dynamic testing of PVB and SG interlayer materials at a range of strain rates and temperatures are presented.For each testing group, five valid tests were successfully completed.Subsequently, average engineering stress-strain curves were produced from the five representative samples.Average parameter values of yield stress, yield strain, failure stress, strain at failure, Young's modulus, and tensile toughness for each testing group were calculated and are presented in Table 2.

PVB Interlayer
The engineering stress-strain responses of PVB at three different dynamic strain rates and three test temperatures are shown in Figure 6.The effect of temperature on the nature of the behavior of PVB is apparent: at low temperature, the dynamic response of PVB is elastoplastic; at ambient temperature, the dynamic response of PVB is viscoelastic and generally bilinear; and at elevated temperature, the dynamic response of PVB is hyperelastic.As the temperature increases, the Young's modulus and tensile toughness decrease, and the material becomes more ductile.There is clear yielding and subsequent strain softening in the dynamic response of PVB at low and ambient temperatures, but not at elevated temperature.This is because the elevated temperature of 60 • C is greater than the glass transition temperature of PVB, which is about 50 • C. At the elevated temperature above the glass transition temperature, the polymer enters a rubbery state where the polymer chains become more mobile, allowing for greater molecular rearrangement and chain slippage.

SG Interlayer
The engineering stress-strain responses of PVB at three different dynamic strain rates and three test temperatures are shown in Figure 7. Unlike PVB, the temperature effect did not change the nature of the dynamic response of SG.At all temperatures, the response was generally elasto-plastic, with a very stiff initial response followed by a clear yielding region and a subsequent strain-softening region.However, there was strain-hardening behavior in the large-strain response of SG at the elevated temperature.Temperature had significant impact on the failure strain of the dynamic response of SG.At low temperatures, SG became brittle and had significantly decreased failure strain.As the temperature increased, the material became much more ductile, resulting in failure strain roughly seven times greater at the elevated temperature compared to the low temperature.At the elevated temperature, the strain-softening region following yielding was much less pronounced.
Effects of the strain rate on the dynamic response of SG could be observed at the low and elevated temperatures.As the strain rate increased, the Young's modulus and yield Strain rate effects on the dynamic responses of PVB are also significant.Generally, as the strain rate increases, the yield stress, failure stress, and tensile toughness increase.These effects become more pronounced as the temperature increases.The strain rate predominantly affects the initial region of the dynamic response.As the strain rate increases, the Young's modulus and yield stress also increase.It was found that strain rate had a negligible effect on the failure strain at low and ambient temperatures; however, the effect was significant at elevated temperatures.At elevated temperatures, the lower speed of the intermediate strain rate of 2 s −1 allowed sufficient time for the molecular chains to rearrange, resulting in more gradual deformation behavior and significantly greater plastic strain than at higher strain rates.

SG Interlayer
The engineering stress-strain responses of PVB at three different dynamic strain rates and three test temperatures are shown in Figure 7. Unlike PVB, the temperature effect did not change the nature of the dynamic response of SG.At all temperatures, the response was generally elasto-plastic, with a very stiff initial response followed by a clear yielding region and a subsequent strain-softening region.However, there was strain-hardening behavior in the large-strain response of SG at the elevated temperature.Temperature had significant impact on the failure strain of the dynamic response of SG.At low temperatures, SG became brittle and had significantly decreased failure strain.As the temperature increased, the material became much more ductile, resulting in failure strain roughly seven times greater at the elevated temperature compared to the low temperature.At the elevated temperature, the strain-softening region following yielding was much less pronounced.

Numerical Modeling
The three-network viscoplastic (TNV) model is a general-purpose, viscoplastic constitutive model that is capable of capturing the thermo-viscoelastic response of many thermoplastics.In this section, the theoretical formulation of the TNV model is outlined.The methodology for calibrating the TNV model using PolymerFEM [54] software is detailed.Material model calibrations for the thermo-viscoelastic properties of the dynamic responses of PVB and SG interlayer materials are presented.A finite element model for simulating the dynamic tensile tests was developed and used to validate the calibrated mod- Effects of the strain rate on the dynamic response of SG could be observed at the low and elevated temperatures.As the strain rate increased, the Young's modulus and yield stress increased, resulting in greater tensile toughness.The strain rate had a greater effect on the Young's modulus of SG at the low temperature than at the elevated temperature.

Numerical Modeling
The three-network viscoplastic (TNV) model is a general-purpose, viscoplastic constitutive model that is capable of capturing the thermo-viscoelastic response of many thermoplastics.In this section, the theoretical formulation of the TNV model is outlined.The methodology for calibrating the TNV model using PolymerFEM [54] software is detailed.Material model calibrations for the thermo-viscoelastic properties of the dynamic responses of PVB and SG interlayer materials are presented.A finite element model for simulating the dynamic tensile tests was developed and used to validate the calibrated models at different temperatures and strain rates.

Material Constitutive Models
The calibrated TNV model for the material comprised three parallel networks, as depicted in Figure 8 in rheological terms.The total deformation gradient F affected all three networks equally, F = F A = F B = F C , and the total Cauchy stress σ was the sum of stresses from each network σ = σ A + σ B + σ C .To simplify the equations in this section, we have omitted the subscripts that differentiate between the networks.This simplification is valid because networks A, B, and C all follow the same constitutive functions.The viscoplastic deformation gradient acting on each network can be separated into elastic and viscoplastic components [54].

2024, 16, x FOR PEER REVIEW
The calibrated TNV model for the material comprised thr picted in Figure 8 in rheological terms.The total deformation networks equally, F = FA = FB = FC, and the total Cauchy stres from each network σ = σA + σB + σC.To simplify the equatio omitted the subscripts that differentiate between the networks because networks A, B, and C all follow the same constitutive deformation gradient acting on each network can be separated components [54].The Cauchy stress acting on a network selected is given by the Yeoh model [53], with Cauchy stress [54]: The Mullins damage factor η A is expressed as follows [54], where σ represents the Cauchy stress, f ε denotes the plastic flow evolution damage factors, [C 10 , C 20 , C 30 ] are parameters of the Yeoh model, J e stands for the determinant of the elastic deformation gradient, η represents the Mullins damage factor, b * represents the left Cauchy-Green deformation tensor, I * 1 denotes the first principal invariant of the Cauchy tensor, and [κ 1 , κ 2 , κ 3 ] are the bulk moduli.Distortional scalars and tensors are indicated with a superscript (*).
where [r, Û, β] are material parameters.η quantifies the softening effect observed in elastomers and other crosslinked rubber-like materials during the initial loading cycles.The Mullins effect, also known as Mullins damage, occurs when an elastomer becomes slightly softer after the first few load cycles, with the mechanical response becoming more repeatable in subsequent cycles.The factor η is used in constitutive models to describe the extent of this softening.The variable U A dev is the current deviatoric strain energy density in network A, and U maxA dev is the max deviatoric strain energy density that network A has been exposed to.The Mullins damage calculation is deactivated if r = 0.The plastic flow evolution damage factor f εA is given by the deferential equation [54]: Here, [ f , c e ] are material parameters, f εA = 1 at t = 0, and .γ v A is a constant introduced for dimensional consistency.The velocity gradient of the viscoplastic flow of network A can be written as [54]: The effective stress driving the viscoplastic flow is defined by the Hill deviatoric stress [54], where [ τA , m] represent material parameters.
. γ 0 = 1/s is a constant introduced for dimensional consistency.R(x) = (x+|x|)/2 denotes the ramp function; τA stands for the flow resistance of the network; mmm represents the stress exponent; f p indicates the pressure-dependent viscoplastic flow; N A specifies the direction of the driving deviatoric stress; and bsign(tr[σ])1 contributes to the volumetric viscoplastic flow, with b being a material parameter.
If the parameter F = 0, then the model yields isotropically, and the effective stress is instead given by the Frobenius norm of the deviatoric stress: In these equations, σ ′ = dev[σ] represents the deviatoric stress.If F < 1 and the other Hill parameters have their default values, then the yield stress will be lower in the 1-direction.If G < 1 and the other Hill parameters have their default values, then the yield stress will be lower in the 2-direction.And similarly, if H < 1 and the other Hill parameters have their default values, then the yield stress will be lower in the 3-direction.The tensor N A specifies the direction of the driving deviatoric stress of the relaxed configuration converted to the current configuration [54].
The pressure dependence of the viscoplastic flow is captured using the function [54]: where p = − 1 3 tr[σ A + σ B + σ C ] is the hydrostatic pressure.In summary, if p 0 > 0, then the yield stress is reduced in a tensile pressure state and increased in a compressive pressure state, and if p 0 < 0, then the yield stress is only reduced in a tensile pressure state.The yield evolution function f εp specifies how the flow resistance evolves with the plastic strain.That dependence is captured using the following equation [54]: where ε p is the Mises strain obtained from .F v A , and [ f f , ε f ] are material parameters.All the TNV model parameters were calibrated based on the experimental data using the PolymerFEM software MCalibration 7.3.0.The calibration procedure flowchart, including the finite element simulations, is illustrated in Figure 9.To begin, the smoothed stress-strain responses from uniaxial tensile tests are input into MCalibration 7.3.0as load cases.When utilizing hyperelastic models, it is critical to select the appropriate model based on the type of experimental data.Hyperelastic models may rely on the first stress invariant, I 1 , such as the neo-Hookean, Yeoh, and Arruda-Boyce models, or not, as in the Moony-Rivlin or Ogden models.This distinction in model classes determines the required dataset.Specifically, non-I 1 -based models necessitate experimental tests across multiple loading modes, while a single loading mode may suffice to calibrate I 1 -based models.Any previously calibrated or known material properties, such as Poisson's ratio or Young's modulus, can also be specified as load cases.
Next, the network types within the TNV model are selected from the PolyUMod library database depending on the type of polymer.Up to three different network types can be employed.The model supports the following hyperelastic components: Yeoh, anisotropic HGOB hyperelastic HGOB (Holzapfel-Gasser-Ogden-Bergstrom), and the hyperfoam model.Each component can be purely elastic or augmented with isotropic or anisotropic power flow and/or the Mullins damage criterion.Hyperelastic component selection is guided by material category: either thermoplastics, elastomers, thermoplastic elastomers (TPEs), or polymer foams.based on the type of experimental data.Hyperelastic models may rely on the first stress invariant, I1, such as the neo-Hookean, Yeoh, and Arruda-Boyce models, or not, as in the Moony-Rivlin or Ogden models.This distinction in model classes determines the required dataset.Specifically, non-I1-based models necessitate experimental tests across multiple loading modes, while a single loading mode may suffice to calibrate I1-based models.Any previously calibrated or known material properties, such as Poisson's ratio or Young's modulus, can also be specified as load cases.Next, the network types within the TNV model are selected from the PolyUMod library database depending on the type of polymer.Up to three different network types can be employed.The model supports the following hyperelastic components: Yeoh, anisotropic HGOB hyperelastic HGOB (Holzapfel-Gasser-Ogden-Bergstrom), and the The calibrations of the PVB and SG models were performed using a two-network configuration.The third network of the TNV model was deselected.The standard Yeoh hyperelastic model was chosen for Network 1 to capture the large strain responses of the materials.After the initial linear regions of the dynamic responses of PVB and SG, the general large strain behavior could be well defined by a simple hyperelastic model.However, viscoplastic thermoplastics present a couple of notable challenges: (1) the Young's modulus is highly dependent on the strain rate; and (2) the large strain responses at different strain rates sometimes converge, like for the dynamic responses of PVB at ambient temperatures (Figure 10b), and at other times, the large strain responses at different strain rates do not converge, as in the case of SG at an elevated temperature (Figure 11c).Therefore, the Yeoh + Power Law Flow Cessation (PSC) model was selected for Network 2. In this model, the power exponent m is modified to be strain rate-dependent.The typical flow rule is given by the equation: .
Polymers 2024, 16,1870 this model, the power exponent  is modified to be strain rate-dependent.The typical flow rule is given by the equation:    The modified flow rule with a strain-dependent -parameter is given by: where the effective power exponent  is defined by: The modified flow rule with a strain-dependent m-parameter is given by: .
where the effective power exponent m e f f is defined by: In this formulation, the power exponent m e f f is controlled by two parameters, m i and m f , which can be optimized to fit the level of strain rate dependence of the Young's modulus.The flow rate is further optimized to control the large strain behavior by inclusion of the flow cessation factor f v .In this approach, the flow rate is given by the following equation: .
γ 0 = 1/s is a constant introduced for dimensional consistency, and the flow cessation factor f v is given by: The β parameter specifies the normalized molecular orientation angle.Note that, with this definition, if A = 0, then the flow cessation is deactivated, and if A = 1, then the flow cessation becomes the same as the equation proposed by Dupaix and Boyce.In other words, A can be used to scale how strong the flow cessation should be.
After the selection of network types, a couple of intrinsic parameters needed to be input.The mass densities of PVB and SG were set to 0.95 and 1.07 g/cm 3 , respectively [24].The Poisson's ratios for PVB and SG were set as load cases to 0.495 and 0.45, respectively [25,26].
Optimization can be carried out using several different numerical methods, such as a random search, a genetic search, a global optimum search, etc.After selection of the material model and network types, initial parameter predictions are generated using the MCalibration 7.3.0software.The parameters which need optimization are selected by the user, and lower and upper bounds can be modified or deselected for each parameter.Once the initial parameters were selected and bounds were set, an automatic extensive optimization was performed, which incorporated an initial random search, the Levenberg-Marquardt search, and the new unconstrained optimization with quadratic approximation (NEWUOA) algorithm.The optimizations performed in this study targeted the normalized mean absolute difference, although it is important to note that several statistical metrices, including the mean square difference or the coefficient of determination, can be utilized.For SG, a fitness weight function of f (x) = 5exp(−x) + 1 was utilized to prioritize the fit of the initial response and softening region.The calibrated models were then exported into Ansys LS-DYNA [62] format and applied to the geometry as user-defined material models.
The calibrations for PVB and SG at each temperature are shown below in Figures 10  and 11, respectively.For each calibration, the normalized mean absolute difference (NMAD) fitness function was chosen.The calibrations resulted in good fits to the experimental data, with the average NMAD errors for PVB and SG being 6.6% and 3.1%, respectively.The calibrated model parameters for PVB at 0 • C, 23

Finite Element Model Components
A finite element computational model was built in the Ansys LS-DYNA to simulate the tensile test under displacement control conditions.The tensile test specimen geometry and boundary conditions are shown in Figure 12a.All nodes of the top grip shown in yellow were restrained in the transverse directions and permitted to move only in the longitudinal direction.All nodes in the bottom grip shown in blue were fixed in all directions.The displacement was applied to the top grip, which caused the specimen to deform in the vertical direction only.Hexahedral eight-node linear elements (HEX8) with 0.635 mm mesh sizes were used (Figure 12b).A convergence study was conducted, and the model converged at a 0.635 mm element size.The PolymerFEM software was integrated into Ansys LS-DYNA to define the material models and input all calibrated material model parameters.

Material Model Validations and Discussion
In order to validate that the calibrated three-network viscoplastic (TNV) models were capable of simulating the thermo-viscoelastic responses of PVB and SG, finite element simulations of tensile tests at the experimental test temperatures and strain rates were performed.The simulations were performed using displacement control and isothermal temperature fields.To accurately replicate the experimental tests, two nodal displacement probes were used to measure the longitudinal displacement of the central 25.4-mm (1-inch) gage length, and the engineering strain history was derived from these nodal displacements.The dynamic force reaction was measured at the bottom fixed support, which is analogous to the load cell used in experimental testing, and the engineering stress history was calculated as the force reaction divided by the original cross-section of the central gage length.For each temperature, simulations were performed at strain rates of 2 and 45 s −1 .
The comparisons between the finite element analysis results and the experimental testing results for PVB and SG are shown below in Figure 13 and Figure 14, respectively.The models predictions agreed well overall with the experimental data and were able to capture substantial changes in behavior due to temperature and strain rate.The fitness of each prediction was characterized by the normalized mean average difference (NMAD) percentage between the finite element analysis test results and experimental data, given in Table 5.The NMAD for each test group was less than 10%.The NMAD was highest for the PVB at 60 °C group, which was similar to the calibration, due to the extreme strain rate effect on the response of PVB at temperature greater than the glass transition temperature.The material model for SG accurately captured the post-peak softening behavior at 0 and 23 °C and the subsequent hardening response of SG at 60 °C.Overall, the constitutive model predictions demonstrate the suitability of the calibrated material models to simulate the thermo-viscoelastic dynamic responses of PVB and SG.

Material Model Validations and Discussion
In order to validate that the calibrated three-network viscoplastic (TNV) models were capable of simulating the thermo-viscoelastic responses of PVB and SG, finite element simulations of tensile tests at the experimental test temperatures and strain rates were performed.The simulations were performed using displacement control and isothermal temperature fields.To accurately replicate the experimental tests, two nodal displacement probes were used to measure the longitudinal displacement of the central 25.4-mm (1-inch) gage length, and the engineering strain history was derived from these nodal displacements.The dynamic force reaction was measured at the bottom fixed support, which is analogous to the load cell used in experimental testing, and the engineering stress history was calculated as the force reaction divided by the original cross-section of the central gage length.For each temperature, simulations were performed at strain rates of 2 and 45 s −1 .
The comparisons between the finite element analysis results and the experimental testing results for PVB and SG are shown below in Figures 13 and 14, respectively.The models predictions agreed well overall with the experimental data and were able to capture substantial changes in behavior due to temperature and strain rate.The fitness of each prediction was characterized by the normalized mean average difference (NMAD) percentage between the finite element analysis test results and experimental data, given in Table 5.The NMAD for each test group was less than 10%.The NMAD was highest for the PVB at 60 • C group, which was similar to the calibration, due to the extreme strain rate effect on the response of PVB at temperature greater than the glass transition temperature.The material model for SG accurately captured the post-peak softening behavior at 0 and 23 • C and the subsequent hardening response of SG at 60 • C. Overall, the constitutive model predictions demonstrate the suitability of the calibrated material models to simulate the thermo-viscoelastic dynamic responses of PVB and SG.

Conclusions and Future Perspectives
In this study, the thermo-viscoelastic dynamic behaviors of PVB-and SG-laminated glass interlayer materials under monotonic loading conditions were investigated experimentally.A drop weight testing machine was modified to perform dynamic testing of the polymers at high strain rates of 20 and 45 s −1 .A pneumatic testing machine was fabricated to perform intermediate strain rate testing at 2 s −1 .Both testing machines were outfitted with a temperature enclosure to perform isothermal testing at temperatures of 0, 23, and 60 • C. The experimental engineering stress-strain data were used to calibrate temperatureand strain rate-dependent material models for PVB and SG using the three-network viscoplastic (TNV) constitutive model.A representative finite element model of the tensile tests performed was constructed, and simulations were performed at each temperature and multiple strain rates to validate the calibrated material models.The main findings of this study are as follows:

•
Strain rate predominantly affected the initial responses of PVB and SG.Initial modulus and yield stress increased with strain rate, resulting in increased toughness.• Temperature affected the characteristic nature of the response of PVB, but not for SG.As temperature increased, the response of PVB transitioned from elastoplastic to bilinear to hyperelastic.At all temperatures, SG interlayer specimens exhibited a linear response followed by yield; strain softening; plastic flow; and, finally, strain hardening at large strains, achieved at elevated temperatures.

•
In general, PVB and SG became stiffer and less ductile as the temperature decreases.The temperature effects were more pronounced in SG than in PVB.For a strain rate of 45 s −1 , as the temperature decreased from 23 • C to 0 • C, the failure strain of SG decreased by 69% and the initial modulus increased by 418%, whereas for PVB, the failure strain decreased by 26% and the initial modulus increased by 200%.• There was clear yielding and subsequent strain softening in the dynamic response of PVB at 0 • C and 23 • C, but not at 60 • C. The elevated temperature of 60 • C was greater than the glass transition temperature of PVB of about 50 • C, and therefore, the behavior was dominated by the viscous component of the viscoelastic response.

•
Temperature-and strain rate-dependent material models for the dynamic responses of PVB and SG interlayer materials were calibrated using the three-network viscoplastic (TNV) model, with average normalized mean absolute differences (NMADs) between the experimental data and model predictions of 6.6% and 2.7% for PVB and SG, respectively.• The FEA results based on the calibrated material models were successfully able to predict the thermal viscoelastic dynamic behaviors of PVB and SG interlayer materials under tensile loading, with average NMADs between the FEA results and experimental data of 6.0% and 4.9% for PVB and SG, respectively.
In the future, additional investigations are needed to improve the material models for their particular use of modeling laminated glass systems.The calibrations at each temperature can be used to develop combined multi-temperature material models to predict the response of the material at any temperature or strain rate within the bounds of the experimental data.Additionally, future work will focus on incorporating an analytical approach similar to the WLF model to account for the temperature dependence of material properties.This enhancement aims to provide a more comprehensive understanding and accurate prediction of the material behavior under varying thermal conditions.During the initial response of laminated glass before glass cracking, the interlayer material acts to transfer shear force between the glass panels during bending.Therefore, the addition of simple shear test data at multiple temperatures to the calibrations can improve the ability of the material models to predict the initial bending response of laminated glass.

Figure 3 .
Figure 3. Pneumatic test setup: (a) pneumatic testing machine with temperature enclosure, select frames from video of an intermediate strain rate test on PVB: start of testing (left) a test (right).

Figure 3 .
Figure 3. Pneumatic test setup: (a) pneumatic testing machine with temperature enclosure, and (b) select frames from video of an intermediate strain rate test on PVB: start of testing (left) and mid-test (right).

Figure 5 .
Figure 5. Temperature control system with drop weight testing machine.

Figure 5 .
Figure 5. Temperature control system with drop weight testing machine.

Figure 5 .
Figure 5. Temperature control system with drop weight testing machine.

Polymers 2024 ,
16, x FOR PEER REVIEW 11 of 26stress increased, resulting in greater tensile toughness.The strain rate had a greater effect on the Young's modulus of SG at the low temperature than at the elevated temperature.

Figure 8 .
Figure 8. Rheological representation of the constitutive three-networ
• C, and 60 • C are given below in Table 3.The calibrated model parameters for SG at 0 • C, 23 • C, and 60 • C are given below in Table 4.

Table 5 .
Fitness of the calibrated TNV material models for PVB and SG5000 given by normalized mean absolute difference (NMAD) (%) for each testing group.